Mechanism of pseudogap probed by a local impurity 
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The response to a local strong non-magnetic impurity in the pseudogap phase is examined in two 
distinctly different scenarios: phase-fluctuation (PF) of pairing field and d- density-wave (DDW) 
order. In the PF scenario, the resonance state is generally double-peaked near the Fermi level, 
and is abruptly broadened by vortex fluctuations slightly above the transition temperature. In the 
DDW scenario, the resonance is single-peaked and remains sharp up to gradual intrinsic thermal 
broadening, and the resonance energy is analytically determined to be at minus of the chemical 
potential. 
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Aside from others, ||l[ two distinctly different scenarios 
are proposed for the pseudogap in cuprates, |^ depend- 
ing on whether the pseudogap phase is independent of the 
pairing gap. In the phase-fluctuation (PF) scenario, ||] 
it is speculated that the normal state contains preformed 
Cooper pairs, and the phase-fluctuation of the pairing 
field destroys superconductivity. As the pairing gap has 
a d-wave symmetry in the internal momentum space, the 
d- wave-like dispersion of the pseudogap Q follows imme- 
diately. An advantage of this scenario is that it involves 
no symmetry breaking, and is adiabatically connected to 
the paramagnetic Mott insulator. Such a normal state is 
not a Fermi liquid. In the second scenario, the normal 
state is free of pairing instability, but is in a symmetry- 
breaking d-density-wave (DDW) state. (| The latter is 
an ordered state of staggered orbital current, and was 
discussed in other contexts already in the early stage of 
high-Tc physics. It creates four hole-like Fermi pock- 
ets in the nodal directions. The volume enclosed by the 
Fermi pockets scales exactly as the doping level x. Thus 
the pseudogap is from the band structure effect. The 
normal state is a Fermi liquid, namely, a DDW metal. 

In this Letter we discuss the resonance state due to 
a strong nonmagnetic local impurity in the pseudogap 
phase, which turns out to be markedly different in these 
scenarios. In the PF scenario, it is two-peaked near the 
Fermi level, broadened by the vortex fluctuations (in ad- 
dition to the intrinsic thermal broadening), and is thus 
strongly temperature dependent near the superconduct- 
ing transition temperature Tc- (This is a complemen- 
tary result to the extended impurity one. In con- 
trast, we analytically verify that the resonance state is 
single-peaked, remaining sharp and pinned at minus of 
the chemical potential in the DDW metal, jsj and iden- 
tify the underlying mechanism for the pinning effect. We 
propose to measure the temperature and doping depen- 
dence of the impurity state by, e.g., scanning tunnelling 
microscopy, |^ in order to tell the mechanism of the pseu- 
dogap (or whether the normal state is a Fermi liquid). 



In a d-wave superconductor, a local strong impurity 
is known to give rise to a resonant state near the Fermi 
level. H^] The state is almost real as the scattering rate 
into the continuum is limited by the vanishing density 
of states near the Fermi level (in an unperturbed sys- 
tem) because of the d-wave pairing symmetry. Qualita- 
tively similar resonant states are found numerically in a 
d-wave superconductor with DDW order. Q This is be- 
cause 1) The resonant impurity state near the Fermi level 
is a generic feature of d-wave pairing, and 2) The only 
significant effect of DDW is to generate a specific band 
structure, on top of which d-wave pairing occurs. 

Returning to the normal state, it is natural to expect 
qualitatively different responses to local impurities in the 
PF and DDW scenarios, since in the DDW normal state 
no pairing occurs. An earlier attempt to address the res- 
onance in the pseudogap phase was made in Ref. |pl]] , but 
with only limited success for PF. Recently the extended 
impurity was discussed in the PF scenario, 0| and a nu- 
merical study was performed for a local impurity in the 
DDW scenario. In this Letter we compare the behav- 
iors of the same local impurity in both PF and DDW 
scenarios. 

Phase- fluctuation scenario: The effective mean field 
hamiltonian in a square lattice for a d-wave supercon- 
ductor may be written as if = + ^■'^■) ~ 

ME^*lc^3*^, where = {lih f^V is the Nambu 
spinor, /i is the chemical potential, 1T3 is the third Pauli 
matrix, and 

with Aij = AoTjij exp{iipij), where Aq is the pairing am- 
plitude, rjij = 1 (—1) for x-dircction (y-dircction) bonds, 
and ipij is the phase. In the PF scenario, the pairing field 
is disordered by thermal and/or quantum fluctuations of 
vortices at zero applied magnetic field. In the follow- 
ing discussion we assume that the vortex fluctuations are 
thermal. 
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It is possible to make a singular gauge transform 
"iii Q-^'Pi'^if^i^i^ so that Aij Ao?7y no longer carries 
the phase, whose effect migrates to the hopping part in 
h.j-. -tas -tase'^'^'-'*'^^'^^/^. In the continuum limit, 
the phase difference between neighboring sites translates 
to the phase gradient 2qs = e~*'^^e*"*, and corresponds 
to the superfluid velocity. It varies at the length scale of 
the London penetration depth, being much larger than 
the Fermi wave-length. Thus it is valid to adopt the 
quasi-classical approximation, in which Fermions see a 
microscopically constant qg while q^ itself varies macro- 
scopically. In this sense, the (matrix) Greens function Gq 
for the ^E" Fermions influenced by q^ is determined by, in 
momentum space, 

^ {iUn - qs • Vk)cro - EkCTS + AkCTl (1) 

Here — —2t{coskx + cos ky) — fi, Ak — 2Ao(cosfc2; — 
cosfcy), Vk = Vk£k- As is well known, the low energy 
excitations in this system are located around the four 
nodes k„=i.2,3.4 = {±K,±K) in the momentum space, 
with —At cos K = fi. The second line in Eq.(^ is the 
usual Doppler approximation. [|l2| At low energies, the 
Doppler shift q^ • Vk can be well approximated by its 
value at the four nodes Z)„ = q^ • Vk„ . We shall use the 
first (second) line of Eq. (|l|) for numerical (qualitative and 
analytical) discussion. [|l3[ 

The real space Greens function is obtained 
by the Fourier transform, G'o(j, j, q^) = 
J2k Go(k, iw„; q) exp[ik • (r^ — r^)]. Of special interest is 
the local Greens function g(iw„; q^) = 6*0(1 = j, zcj„; q^). 
Due to the d-wave pairing symmetry, 912 = 521 — 0. 
(We suppress the arguments if apphcable.) On the other 
hand, 

gn{iuJn;qs) = -^22 = J dENo{E;(is)/{iuJn - E), (2) 

4 

No{E; q,) ~ (1/8) ^ 51 - '^Dr.l/i^^tAo) 

n— 1 z^— it 

- ^i{E - iyD„)hgniE - vDr,)/K% (3) 

with A = A[Kt^All{e + Ag)]i/4. A cutoff at \E\ > E^^ 
min(4<,4Ao) is necessary in applying Eq.(^. Anticipat- 
ing its effect in the impurity scattering, we point out 
briefly the behavior of A^o^ 1) It exhibits a four- fold 
symmetry in the direction of q^. 2) At \E\ ^ Ec and 
\Dn\ <C Ec, the leading contribution comes from the first 
term in Nq. It is particle-hole symmetric around E = Q 
at any and qs. This is the fundamental property of 
d-wave pairing between time-reversed particles. 3) Away 
from half filling (/x ^ 0), there is a slight asymmetry. To 
the first order in /x, this is included in the second term in 
Nq. Clearly, Nq{E — 0;qs) oc q^, as first pointed out by 

Voiovik. |T|i 



In the presence of a local scattering potential at site 
i = 0, the corresponding Greens function G can be ob- 
tained within the T-matrix approximation (which is ex- 
act if the impurity does not spoil the pairing field), 

G{i,3,iu)n;<=^s) = GQ{i,3,iujn](\s) 

+Go(i, 0, icj„; (\s)T{iujn] qs)Go(0, j, zw„; q^), (4) 
T~'^{iujn] qs) = l/iVa^ + VmCQ) - g{iuJn; qs), (5) 

where V [Vm) is the nonmagnetic (magnetic) scattering 
strength. To simplify our discussion, we shall consider 
nonmagnetic scattering only. With the impurity, the lo- 
cal density of states (LDOS) is site dependent. At site i 
it is given by 

N{i,uj;cis) = -(l/7r)ImGii(i,i,u; + fO+;qs). (6) 

We emphasize that the off-diagonal elements of Go in 
Eq.(^ contribute to the density of states. Conceptually, 
neglecting such contributions, as in Ref. |jri|, the theory 
would be hardly related to pairing. 

Let us discuss the qualitative behavior of the LDOS in 
our case. Since g is diagonal, so is the T-matrix. After 
analytical continuation iujn — > -I- iO^, the resonance in 
LDOS is determined by Re(rfi^) = or Re(T2^^) = 0. 
This is equivalent to Ke[gii{±uj + iO'*";qs)] = V~^, and 
immediately implies two resonance peaks in a general 
situation, an essential feature of pairing. For a strong 
scatter, — > 0. In the case of = 0, a sharp res- 
onance exists in the LDOS at the sites nearest to the 
impurity. [|o| The behavior at ^ is as follows: 1) 
For jjL — {) and = 0, we have perfect particle-hole 
symmetry, so that Re[(7ii(iO+; q^)] = from Eq.(||). The 
resonance is at a; = 0. However, it is not sharp, and 
its width scales with iV(0;qs) (oc q^) in accordance with 
the Fermi golden rule. 2) The effect of a finite ^ and/or 
is to generate a slight asymmetry, and thus splits 
the resonance into two peaks, situated on either side of 
the Fermi level. Their widths are identical (different) if 
II = {fi 0) because of the behavior of A'o. Moreover, 
at large enough Qs the splitting will be smeared due to 
the broadening of both peaks. 

Before we proceed, we predict from the above results 
that even below Tc the resonance may be broadened by 
an in-plane transport current and/or a nearby static vor- 
tex. Since A'o is four-fold symmetric in the direction 
of the relevant super-current, so is the broadening phe- 
nomenon. 

Thermal phase fluctuations are governed by the 
Kosterlitz-Thouless (KT) theory. Q In the quasi- 
classical approximation, this just amounts to an average 
over q.,. For illustrative purpose, one can assume a Gaus- 
sian distribution exp(— qg/2n^)/(27rn^) for q.,, with n„ 
scaling with the density of free vortices (anti- vortices). 
The average LDOS is thus N{i,uj) = (Af(i, w; q^)). 
The averaging makes it inconvenient to determine the 
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resonance exactly, since it is the Greens function that 
should be averaged instead of the T-matrix alone. How- 
ever, qualitative features of the resonance in the LDOS 



are roughly the same as discussed above, but with a char- 
acteristic scale of Qs given by y/n^, the inverse vortex 
spacing. 
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Results with Ao = OATt, 



FIG. 1. Results with Ao = OATt, fi = -0.3t and V = lOOt. (a) iV(r„„,a;) versus u. Solid lines: n„ = ~ 10"*^, 10""*, 10"'', 
and 5 x 10"'^ with decreasing peaks. The dotted line is the LDOS at rii, = and V = for comparison, (b) A^(r, 0.05t) at 
0. The impurity is at the center, (c) The same as (b) for riv — 5 x 10~^. The gray scale is the same in (b) and (c). 



All of the above predictions are indeed seen in our nu- 
merical results. We use Eqs.(0) (the first line) and (^- 
and average over q^. In Figs.l we present results 
for Ao = O.m, n = -OM and V = lOOt. (The re- 
sults do not change much in the unitary limit V :s> t.) 
Fig. 1(a) plots the LDOS at r„„ nearest to the impu- 
rity as a function of energy. For < 10"^ (or vor- 
tex spacing > 10^ in units of crystal lattice constant 
ao), the resonance is sharp and indistinguishable from 
that with no vortices at all. However, it begins to de- 
grade at Tiv > 10""* {dy < lO^flo), and becomes almost 
featureless at n„ = 5 x 10"'^ {d^ = 14ao). Translat- 
ing to the temperature dependence using the KT expres- 
sion n„ ~ exp[— aTc/ (T — Tg)] (with a = 5 for estima- 
tion), one expects no significant change of the resonance 
at T—Tc < 0.07Tc, but it is suddenly degraded as soon as 
T-Tc > O.lTc. For = 40K as in a typical underdoped 
cuprate, the temperature window for this phenomenon to 
happen is within AK. While the exact number should not 
be taken seriously, the sudden degrading of the resonant 
impurity state is a robust and peculiar feature of the PF 
scenario. A similar case was found for an extended im- 
purity elsewhere. The dotted line in Fig. 1(a) is the 
DOS at rit, = and V = ioi comparison. It also 
shows the slight doping-induced particle-hole asymme- 
try. Fig. 1(b) shows the spatial distribution of LDOS at 
uj — 0.05i (<C |/i|), one of the resonance energies, when 
ri„ = 0, |l^ which should be compared to Fig. 1(c) at 
n„ = 5 X 10"'^, upon which the contrast for the four-fold 
structure is much weaker. 



DDW scenario: We assume that the effective Hamil- 
tonian for the uniform DDW metal is |^ 



H 



it- 



iDi-j)clc. + h.c. 



(7) 



Doi]i.j{-iy with 



where spin index is suppressed, Dij 
Dci being the DDW order parameter. 

It proves useful to introduce two sublattices A and B, 
and denote — A-i and Cjgs = Bi. The real-space 
sublattice Greens function is 



^BA ^'^ 



'BB 



E 



Ak,t.exp(zk ■ Tij) 



(8) 



where Ak., 



\/X^ + D^, = 2t{coskx + cosfcy), and Z3k = 

2Dq{cos kx — cos ky). Note that the summation over the 
momentum space is limited to the reduced Brillouin zone. 
In order to study the impurity problem, wc need the real- 
space Greens function G^c'\i,j, icun) in terms of the orig- 
inal c-electrons. This is related to the above as 



(9) 



where a,(3 = A, B. 

The unperturbed on-site Greens function is inde- 
pendent of sublattices, gditOn) = Gc""* (0, 0, zw„) = 
/ dENo{E)/{iu;„-E) with NoiE) = ^^.^i 5{E-vEk + 
^i) \E + ^i\/{STTDt) being the unperturbed DOS. The 
second equality in A^o requires a cutoff a.i \E + ^\ > Ec = 
min(4£',4t). The symmetry around E — in A^Oj in- 
stead of at i? = in the case of pairing, is exact in 
our model. This is because doping the system does not 
change the two symmetric bands generated by DDW, but 
just shift the Fermi level. 

Now, the Greens function in the presence of the impu- 
rity can again be obtained within the T-matrix formula- 
tion. 
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Gc{i,j,iujn) = GfHi,j,iuJn) 

+G(") (^, 0, ^c^„)T(^c^„)G(°) (0, j, iun), (10) 
T-\tu;n)^V-^ -gdiiOn). (11) 

The resonance state is defined by Ke[T^^{u! + iO+)] = 0, 
or equivalently Re[gc{uJ + iO+)] — V~^. Using the ap- 
proximate iVo: 3c(^ + iO^) is given by 



-sgn(a; + /i)iVo(w)ln[ 



l]-inNo{Lu). (12) 



Thus the resonance occurs below (above) —fi for a finite 
positive (negative) potential V . For a nonmagnetic im- 
purity, it is single-peaked because of the unique condition 
for the resonance to occur. Again the energy width of the 
resonance scales with Nq. In the unitary limit V^^ 0, 
the resonance energy is w = — /i from the above expres- 
sion of gc and iVo(~M) — 0. In fact this result is exact 
since the exact symmetry in Nq mentioned above guar- 
antees Re[5c(— M + *0+)] = 0. This resonance energy is 
exactly at the mid-point of the two symmetric bands, 
in much the same way as the mid-gap state exists in a 
semiconductor. Furthermore it would be infinitely sharp 
since No{—fJ.) = 0. 

The DDW order should not fluctuate significantly once 
it is well developed, because it is an Ising-like order pa- 
rameter so that no Goldstone mode exists. Therefore 
the thermal rounding of the resonance is gradual, with 
no abrupt change just above Tc, in contrast to the case 
in the PF scenario. 
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FIG. 2. Results with Do = 0.17t and = 
A'"(rnn,'^) versus uj. (b) N{r,uj — —/j.) versus r. 



-0.3f. 



The LDOS N{i,uj) can be easily calculated from 
Eq. (|l0|) . In Fig. 2(a) we present the LDOS at r„„ near- 
est to the impurity. The resonance is single-peaked and 
robust as long as t,D. In Fig. 2(b) we present the 

spatial dependence of the LDOS at the resonance en- 
ergy. The pattern is also four-fold symmetric, similar 
to the case of a d-wave superconductor except that the 
resonance here is at — instead of at zero energy. 

After the submission of this paper, we become aware 
of a related but independent work. 
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